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We analyze two recent reactions of Belle, producing DD and DD* states that have an enhance- 
ment of the invariant DD, DD* mass distribution close to threshold, from the point of view that 
they might be indicative of the existence of a hidden charm scalar and an axial vector meson states 
below DD or DD* thresholds, respectively. We conclude that the data is compatible with the 
existing prediction of a hidden charm scalar meson with mass around 3700 MeV, though other pos- 
sibilities cannot be discarded. The peak seen in the DD* spectrum above threshold is, however, 
unlikely to be due to a threshold enhancement produced by the presence, below threshold, of the 
hidden charm axial vector meson X(3872). 

PACS numbers: 



I. INTRODUCTION 

The observation of new states with open and hidden 
charm has sparked the interest of many experimental and 
theoretical research groups. In order to accommodate the 
newly found states, many theoretical explanations have 
been suggested and there is still controversy about the 
structure of many of these resonances. New and precise 
data will help understand these states better and it is 
important to test the theoretical models against these 
data. In that sense, the recent reactions of Belle produc- 
ing DD and DD* final particles [l| bring relevant infor- 
mation concerning hidden charm states that we want to 
exploit in the present work. 

The reactions e+e" -> J/ipD^ 5W , with different 
charmed meson pairs, have been recently observed by 
the Belle collaboration [l| . We study here the cases with 
DD and DD* + c.c. pairs in the final state. The invari- 
ant mass distributions for D pair production in these two 
cases present an important enhancement close to the two 
meson threshold, which have led to the claim of two new 
resonance states pj. 

In three recent papers 0, [H, 0] it was shown that 
bound states lying close to threshold for meson produc- 
tion can cause a strong enhancement in the cross sec- 
tion above threshold. Moreover, following similar steps 
as in @, 0, H, [H, we have developed a model for gen- 
erating dynamically resonances with open and hidden 
charm quantum numbers in the scalar and axial sectors 
[H, [j| . The model describes a large number of observed 
resonances as dynamically generated states in agreement 
with d, 0, H, [1] • However, some other states predicted 
with narrow widths in @, 0, H, @ appear in [H, |j| as very 
broad resonances, which could explain why such states 
are not being observed. In particular the extension of 



[H, [1[ to the hidden charm sector predicts new meson 
states both scalar and axial vector. 

The predicted hidden charm resonances are one scalar 
and two axial vectors. One of these axials can accommo- 
date the X(3872) [THHH which most probably has pos- 
itive C-parity conjugation but the other axial pole 
has negative C-parity and is nearly degenerated in mass 
with the positive one. The scalar resonance is predicted 
around 3.7 GeV with a very small width and there is so 
far no experimental evidence for it. Since the scalar pole 
lies close to the DD production threshold and the axial 
ones close to the DD* + c.c. threshold, it is possible that 
the processes measured by Belle are dominated by these 
dynamically generated resonances, so that the threshold 
enhancement is a consequence of these poles instead of 
being due to new resonances above threshold. 

We briefly present, in the next section, our model for 
generating dynamically resonances and show how to ap- 
ply it for calculating cross sections for J plus D pair 
production from electron positron annihilation. Section 
III shows our results and in section IV we present our 
conclusions. 



II. FRAMEWORK 

The reaction e + e~ — > J/ipDD can be described by the 
diagram in fig. [1] if one assumes that the DD pair comes 
from a resonance. 
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FIG. 1: Feynman diagram for the process e + e — ► J/tpDD 
Close to threshold the only part of this amplitude 
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which is strongly energy dependent is the X propaga- 
tor and all other parts can be factorized, so that we can 
write 
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if we describe the X resonance as a Breit-Wigner type. 

The cross section would then be given by an integral 
over the phase space of the three particles in the final 
state: 
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Assuming that T depends only on the DD invariant 
mass, one can evaluate from cq. ([2} the differential cross 
section: 
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where s is the center of mass energy of the electron 
positron pair squared and | k | and \~p f \ are given by: 
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(4) 
(5) 



Where A 1 ^ 2 (s, m 2 , M 2 ) is the usual Kallen function. 

Let us briefly see how the approach of [!, H[ generates 
dynamically some resonances. The model starts by writ- 
ing a Lagrangian for the interaction of the meson fields 
belonging to the 15-plet of SU(4). For the pseudoscalar 
field we take 
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and an analog one for the vector meson field. 

For each one of these fields a current (J p and J^) is 
constructed and the Lagrangians are build by connecting 
these two currents: 
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The Lagrangians in (J6j) and {7]) are SU(4) invariant, 
but since this is not an exact symmetry in Nature, one 
breaks it by suppressing the terms in these Lagrangians 
where the underlying interaction is driven by the ex- 
change of a heavy vector meson. A different SU(4) sym- 
metry breaking pattern borrowed from (l3j was also stud- 
ied in 0] in order to estimate the uncertainties of the 
model. 

The tree level amplitudes for meson meson scatter- 
ing are obtained from the Lagrangians for all possible 



processes. All meson pairs forming a system with equal 
quantum numbers span a coupled channel space. All am- 
plitudes, projected in s-wave, for the processes in each 
coupled channel space are collected into a matrix. We 
call this matrix the potential, V. 

The potential is then used to solve the Bethe-Salpeter 
equation that, in an on-shell formalism [3,[lB|, assumes 
an algebraic form: 



T = (L-VG^V 



T 
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(8) 
(9) 



Equation ([5]) is used for calculating the T-matrix in the 
scattering of two pseudoscalars that, in s-wave have the 
quantum numbers of a scalar resonance. Equation © 
is used in the scattering of pseudoscalars with vector 
mesons, takes into account the polarization of the vec- 
tor mesons and the resulting states have quantum num- 
bers of an axial. The matrix G is diagonal, coming from 
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the loop function of the two meson propagators of the 
intermediate states. 

The loop function has its imaginary part fixed in or- 
der to ensure unitarity of the T-matrix. Resonances are 
identified as poles in the second Riemann sheet of the 
T-matrix. 

With this model many known states have been repro- 
duced: a scalar antitriplet with charm quantum num- 
ber that has been identified with the Z)* (2317) and 
Dq (2400) states and two axial antitriplets identified 
with the D s i(2460), £>i(2430), Aa(2536) and £>i(2420). 
Apart from those, some broad sextets are predicted in the 
scalar and axial sectors, though the ones in the scalar 
sector are probably too broad to have any experimen- 
tal relevance. The light scalars and the light axials are 
also reproduced in our model, and finally in the hidden 
charm sector three poles are generated in the T-matrix, 
one scalar around 3.7 GeV and two axials, nearly degen- 
erated in mass, with opposite parity and close to the DD* 
threshold, one of which is identified with the A (3872). 

Since the dynamically generated states are character- 
ized by poles appearing in the unitary T-matrices, the 
dynamics of our approach is incorporated in the e + e~ — > 
J/ipDD process by substituting the Breit-Wigner ampli- 
tude of (Q} by the DD T-matrix calculated from eq. (jHJ). 
For the reaction e + e~ — > J/ipDD* everything is done 
analogously using eq. © for the T-matrix. 

In the following section we explain how we compare 
our results to Belle's data. 



III. RESULTS 

Belle has measured the differential cross section for 
J/tpDD, J/ipDD* and J/ij)D*D* production from elec- 
tron positron collision at center of mass energy -y/s=10.6 
GeV [l[ . We are going to study the first two cases, where 
the hidden charm states generated in the model of [H, [f| 
could be related to. The Belle's measurement produces 
invariant mass distributions for the DD and DD* that 
range from threshold up to 5.0 GeV. Our model is, in 
principle reliable for energies within few hundreds of MeV 
from the thresholds, so we are going to compare numer- 
ically our results with the data up to 4.2 GeV. 

The experiment measures counts per bin. In the case 
of a DD pair, the bins have 50 MeV width, while for 
the DD* pair they have 25 MeV. To compare the shape 
of our theoretical calculation with the experimental data 
we integrate our theoretical curve in bins of the same size 
as the experiment and normalize our results so that the 
total integral of our curve matches the total number of 
events measured in the invariant mass range up to 4.2 
GeV. 

The comparison is made by the standard \ 2 test. The 
value of x 2 divided by the number of degrees of freedom 
is given by: 



TABLE I: Results of Mx and \ 2 f° r different values of cxh- 





Mx(MeV) 


x 1 

d.o.f 


-1.4 


3701.93-i0.08 


0.38 


-1.3 


3718.93-i0.06 


0.36 


-1.2 


3728.12-i0.03 


0.49 


-1.1 


Cusp 


0.66 



X 2 _ 1 \ ^ iXtheo Yexp) 2 /,„n 

d.o.f. ~ (AY exp )i [ } 

where Y is the number of counts in each bin, AY" is the 
experimental uncertainty in each measurement and N is 
the total number of points. 

As described in in the heavy sector the model has 

one free parameter, an which is the subtraction constant 
in the loop for channels with at least one heavy particle. 
In these previous papers this parameter has been fitted 
so that the pole in the C=l, S=l, 1=0 sector matches the 
observed state with these quantum numbers (D s o(2317) 
and D s i(2460) for scalar and axial states, respectively). 
The channels in this sector have always one heavy and 
one light meson, and in principle one could fit a different 
a for channels with two heavy particles. Here we are go- 
ing to present results for different values of a in channels 
with hidden charm (doubly heavy channels). Since we 
are working with the C=0 sector, we have also channels 
involving only light mesons. These have negligible influ- 
ence in the pole position of the hidden charm poles as 
shown in so we leave at constant. 

In table U we show results, for different values of ajj, 
of the pole position of the hidden charm resonance in 
the scalar sector, and the value of \ 2 calculated with the 
data from Belle, with combinatorial background already 
subtracted, for all points below 4.2 GeV in the J /ifjDD 
production. Fig. [2] shows plots of our theoretical his- 
tograms compared with experimental data. Note that 
although we are plotting all points until 5.0 GeV, only 
the ones below 4.2 have been used in the calculation of 
X 2 and in the normalization of the theoretical curves. 

The x 2 values obtained in table [J are all below 1, indi- 
cating a good fit to the data in all curves. This is in part 
due to the large experimental errors, but the clear mes- 
sage is that the presence of a pole below DD threshold 
is enough to produce the observed enhancement of the 
cross section for this reaction in the DD invariant mass 
above threshold. The results of table U and inspection of 
fig. [2] show some preference for values of a#=-1.3, -1.4, 
which would correspond to the hidden charm scalar with 
mass slightly above 3700 MeV. 

For the production of J/ipDD* we use the model for 
generating axial resonances. In this case the resonance 
A in fig. [1] should be identified with the A" (3872) gener- 
ated by our model. Note that our predicted state with 
negative C-parity does not fit here, since this experiment 
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FIG. 2: Theoretical histograms compared with data from [l|] 
for DD invariant mass distribution. 



FIG. 3: Theoretical histograms compared with data from 
for DD* invariant mass distribution. 



TABLE II: Results of M x and x 2 for different values of an- 



OtH 


M x (MeV) 


x 2 

d.o.f 


-1.40 


3865.09-i0.00 


2.84 


-1.35 


3870.07-i0.00 


3.84 


-1.30 


3872.67-i0.00 


5.04 


-1.25 


Cusp 


5.01 



selects a positive C-parity state for the X. Table HT1 shows 
results for Mx and x 2 f° r different values of an- Since 
the state A (3872) is known and has a rather precise mass, 
we have chosen a smaller range to vary the parameter a 
in order to get the mass of the X closer to its experimen- 
tal value. Fig. [3] compares our theoretical results with 
the experimental data from Belle. 

In this case the x 2 obtained is in all cases of the order 
of 3 or higher, clearly indicating a poor fit to the data. 

The peaks seen in the experiment have been fitted with 
Breit-Wigner like resonances in Q. In the experimental 
work of [l[ the two peaks discussed here have been fitted 
in terms of Breit-Wigner amplitudes, suggesting two new 
resonances. In order to make the results obtained here 
more meaningful, we also perform such a fit and compare 
the results. We take the same Breit-Wigner parameters 
suggested in the experimental paper. The scalar reso- 
nance has Mx=3878 MeV and T x =347 MeV and the 
axial one has M x = 3942 MeV and T x = 37 MeV. We 
show the results obtained by fitting a Breit-Wigner form 
from eq. |T|) in T of eq. ([3|) in figs. [4j [5] Additionally we 
calculate \ 2 and find x 2 /d.o./=1.04 for the DD distribu- 
tion and X 2 /d-o.f=im for the DD* distribution. The 
value of x 2 f° r the DD distribution can be improved if 
we take different parameters for the Breit-Wigner reso- 
nance. Taking for the fit Af x =3750 MeV and T x =250 
MeV we obtain a value of x 2 /<io./=0.55, comparable 



with those obtained in our analysis. The value of x 2 f° r 
the DD* distribution is undoubtedly better in the case 
of a Breit-Wigner fit that in our analysis assuming the 
A (3872) resonance as responsible for the peak. 

As a consequence of the discussion, our conclusions 
would be a support for a new resonance around 3940 
MeV as suggested in [l[ , while for the case of the broad 
peak seen in DD, the weak case in favor of a new state 
around 3880 MeV discussed in [l[ is further weakened 
by the analysis done here, showing that the results are 
compatible with the presence of a scalar hidden charm 
state with mass around 3700 MeV. 
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FIG. 4: Histograms calculated with Breit-Wigner resonance 
with mass Mx=3880 MeV compared to data. 
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FIG. 5: Histograms calculated with Breit-Wigner resonance 
with mass Mx=3940 MeV compared to data. 



IV. CONCLUSIONS 

We have applied our model for generating dynamically 
resonances with hidden charm in order to describe differ- 
ential cross sections in reactions observed by Belle where 
D meson pairs were produced together with J/tjj from 
electron positron collisions at a center of mass energy of 
10.6 GeV. 

In the experiment two peaks are observed, one in the 
DD invariant mass distribution around 3880 MeV and 
another in the DD* invariant mass distribution around 
3940 MeV. We studied here the possibility that these 
peaks might be a threshold enhancement caused by res- 



onances below the threshold for D pair production. 

While our theoretical model describes very well the 
peak in the DD invariant mass as caused by a scalar 
resonance which is mainly a DD bound state, the % 2 val- 
ues obtained in describing the peak in the DD* invariant 
mass distribution are too big, indicating that this peak 
is unlikely to be caused by the A" (3872) resonance re- 
produced in our approach as a dynamically generated 
resonance. 

The results obtained for the invariant mass distribu- 
tion calculated with our model in comparison with data 
are compatible with the existence of the predicted hid- 
den charm scalar resonance below DD threshold, but the 
large experimental errors do not discard other possibil- 
ities. Further experiments to search for this predicted 
states should be most welcome. In particular the radia- 
tive decay of this predicted state leading to J/ i/jj studied 
in [l6| would be a good test for it. 
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